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ABSTRACT 

After having reported the detection of X-rays emitted by the peculiar system HD 5980, we 
assess here the origin of this high-energy emission from additional X-ray observations obtained 
with XMM -Newton. This research provides the first detection of apparently periodic X-ray 
emission from hot gas produced by the collision of winds in an evolved massive binary outside 
the Milky Way. It also provides the first X-ray monitoring of a Luminous Blue Variable only 
years after its eruption and shows that the source of the X-rays is not associated with the ejecta. 

Subject headings: stars: individual (HD 5980) - stars: Wolf-Rayet - stars: winds - Magellanic Clouds - 

1. Introduction 

The most massive stars, which are also the most 
luminous (and for most of their lives the hottest) 
stellar objects, have a large impact on their host 
galaxies. For example, throughout their evolution, 
these stars convert simple into more complex ele- 
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merits, and distribute them into space by a com- 
plex process of mass loss. They are thus largely 
responsible for the chemical enrichment of the Uni- 
verse (Massey 2003). Their intense, hot pho- 
ton emission is also able to ionize the surrounding 
medium (leading to the formation of Hli regions) 
and to drive a continuous ejection of matter in the 
form of a stellar wind. Finally, massive stars in- 
ject a large amount of mechanical energy into their 
host galaxies, through e.g. supersonic winds, large 
transient eruptions (like those of Luminous Blue 
Variables, e.g. 77 Carinae), or GRB/supernova ex- 
plosions (Massey 2003). 

In massive early-type binaries, the supersonic 
outflow from one star collides with that of its com- 
panion. This collision provokes a strong heating 
of the shocked wind, leading to the emission of X- 
rays (see e.g. Stevens et al. 1992). Hints of this 
high-energy phenomenon were already found two 
decades ago, when the first X-ray observatories 
detected overluminosities in hot binaries (Pollock 
1987; Chlebowski & Garmany 1991). However, 
since additional X-ray emission can have several 
origins, not restricted to colliding winds (CW), 
indisputable observational evidence of CW X-ray 
emission had to await the advent of a new genera- 
tion of sensitive X-ray facilities (Chandra, XMM- 
Newton). Indeed, only detailed spectroscopic in- 
vestigations and careful monitoring in the X-rays 
can bring to light the actual properties of the hot 
gas produced by the wind-wind' collision. Most 
notably, phase-locked variations in the X-ray do- 
main are produced by varying separation in eccen- 
tric binaries (that changes the intrinsic strength 
of the collision) or varying line-of-sight opacity 
as the stars revolve around each other: see e.g. 
HD 152248 (Sana et al. 2004), HD 93403 (Rauw 
et al. 2002), WR25 (Pollock & Corcoran 2006), 
and 7 2 Vel (Schild et al. 2004). Up to now, all 
studied X-ray colliding wind (XCW) systems be- 
long to our Galaxy - a few XCW candidates have 
been proposed in 30 Doradus, but solely on the ba- 
sis of an X-ray overluminosity (Portegies Zwart et 
al. 2002), which is insufficient to ascertain their 
nature. Studying the CW phenomenon in other 
galaxies can provide an important probe of the 
mass-loss process in different environments with 
different metallicities. 

HD 5980, the most peculiar massive star in the 
Small Magellanic Cloud (SMC), lies on the periph- 


ery of the large cluster NGC 346 associated with 
the giant Hu region N66. It is a multiple sys- 
tem whose main component, star A, underwent 
two LBV-like eruptions in 1993-1994, increasing 
its brightness by up to 3 magnitudes (Jones & 
Sterken 1997). Together with star B, believed 
to be an early Wolf-Rayet star of the nitrogen se- 
quence, it forms a close eclipsing binary system 
whose period is 19.3d, eccentricity is 0.3 and in- 
clination very close to 90° (Sterken & Breysacher 
1997). Using the ephemeris of the latter authors, 
Star A (resp. B) eclipses its companion at </>— 0 
(resp. 0.36), whereas periastron and apastron oc- 
cur at phases 0.09 and 0.59, respectively. A third 
stellar object, star C (probably an early O-type 
star), contaminates the light of the system, but 
it is still unclear if this star is just a line-of-sight 
coincidence or an object gravitationally bound to 
the close AB pair. The presence of an additional 
component (e.g. a close orbiting neutron star) 
has been proposed, based on a 7-hour periodic- 
ity seen in spectral and photometric variations and 
stochastic polarimetry changes, but is still debated 
(see e.g. Villar-Sbaffi et al. 2003). Moffat et al. 
(1998) infer the presence of a wind-wind collision 
region in the AB system on the basis of strong 
emission-line variability. From UV observations, 
Koenigsberger et al. (2000) and Koenigsberger 
(2004) conclude that the orientation of the shock 
cone is such that it wraps around Star B, con- 
sistent with the notion that Star A possesses the 
more powerful of the two winds. 

The detection of X-ray emission from HD 5980 
was reported for the first time by Naze et al. 
(2002). Variations of the X-ray brightness were 
detected, in the short-term range (Naze et al. 
2002) as well as the long-term range (Naze et al. 
2004). However, these observations were not suf- 
ficient to determine the nature of these variations 
(e.g. variations driven by changes in the shocked 
ejecta from the 1994 eruption, or phase-dependent 
changes due to wind-wind collisions in the AB bi- 
nary), and the origin of the X-ray emission there- 
fore remained uncertain. To resolve this, we re- 
quested additional monitoring of HD 5980 with 
XMM-Newton. 

This letter is organized as follows. Section 2 
describes the observations while Sections 3 and 4 
present the results of the dedicated XMM -Newton 
campaign and their interpretation, respectively. 
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Section 5 give our conclusions. 

2. Observations and Data Reduction 

During the past years, HD 5980 was observed 
once with Chandra during lOOks and five times 2 
with XMM-Newton for approximately 20 ks each 
time (see Table 1). For the latter observations, the 
three European Photon Imaging Cameras (EPICs) 
were operated in the standard, full-frame mode 
(except for the first pn dataset) and a medium fil- 
ter was used to reject optical light. Note that the 
observations were carefully planned to sample cru- 
cial phases of the 19.3d orbit (Fig. 1); moreover, 
the last three observations were obtained during 
the same 19.3d orbit of HD 5980. To ensure a 
coherent reduction, we used the Science Analysis 
System (SAS) software, version 7.0, to re-process 
all X-ray data. After the pipeline chains, filters 
recommended by the SAS team were applied in 
order to keep only the best data (see e.g. Naze et 
al. 2006). For consistency, no additional temporal 
filtering was done. 

To each XMM-Newton dataset, we applied the 
SAS source detection algorithm ( edetect-chain ) in 
a region of 150” radius around HD 5980, in order 
to derive the best value of the centroid of the X-ray 
emission. Since HD 5980 appears surrounded by a 
soft, X-ray bright SuperNova Remnant (SNR, see 
Naze et al. 2002), the detection was restricted to 
the 1.5-10. keV range, to minimize the contamina- 
tion from this soft extended source. As a check, we 
compared these results with the output of another 
algorithm (SAS task eregionanalyse ) that simply 
calculates the number of counts in the region con- 
sidered 3 and corrects them only for the Encircled 
Energy Fraction (about 0.8 in our case). Within 
the error bars, both methods agree with each other 
so we only present the results from the latter (see 
Table 1). 


2 A sixth observation taken during Rev. 0970 was strongly 
affected by a flare, rendering the data unusable. 

3 We used a 25” region centered on HD 5980, together with a 
background region of the same area but offset from HD 5980 
by 22s east in RA (or —2040 px in X) and 82” south in DEC 
(or —1640 px in Y). Using an annular background region 
yields the same results, though noisier. 


3. Results 

It must first be noted that two XMM-Newton 
datasets were taken, intentionally, at the same 
phase (<f> ~ 0.36). Although these two observa- 
tions were obtained during XMM-Newton Revs. 
0157 (2000 Oct.) and 1094 (2005 Nov.), i.e. sep- 
arated by 5 years or 97 orbital revolutions of the 
AB system, they present very similar values for the 
count rate of HD 5980 (Fig. 2 and Table 1). This 
means that post-eruption, epoch-dependent vari- 
ations are minimal. Indeed, if a significant part of 
the X-ray emission came from the collision of the 
fast wind with the slower ejecta associated with 
the 1994 eruption, a monotonic decrease in the 
X-ray brightness and a decline in the X-ray tem- 
perature with time would be expected as this in- 
teraction zone expands and cools. 

Second, coherent variations with phase are now 
clearly detected, with X-ray emission steadily in- 
creasing towards <j> = 0.36 (Figs. 2 and 3). At that 
time, the emission also appears slightly harder, 
though we cannot exclude a constant value of the 
HR at the 2 a level. This bright phase corresponds 
to the eclipse of Star A by its companion, i.e. when 
the opening of the shock cone, slightly concave 
with respect to Star B, is directed towards the ob- 
server (Fig. 1). 

Although we cannot compare directly the 
Chandra and XMM-Newton count rates (because 
of background contamination, different spatial res- 
olution, and cross-calibration problems), it must 
be noted that the variations measured by Chandra 
confirm the trend seen by XMM-Newton. During 
the 100 ks of the Chandra observation, the count 
rate clearly increased from 2 to 4xl0” 3 ctss~ 1 
in the 0.3-10. keV band (Naze et al. 2002) 
or from 1 to 2xl0~ 3 ctss'~ 1 in the 1.5-10. keV 
band (Fig. 2). These data covered phases ranging 
from 0.24 to 0.30. In this interval, the XMM- 
Newton light curve indeed predicts an increase of 
~ 1 x 10““ 3 ctss”* 1 in the count rate. 

Finally, we investigated the X-ray spectrum of 
HD 5980. Of course, the faintness of this dis- 
tant source and the contamination at low ener- 
gies by the superposed SNR prohibit any detailed 
study. However, since the spectral properties of 
the SNR contamination can be determined thanks 
to the high-resolution Chandra data of Naze et al. 
(2002), a first, general spectral evaluation can be 
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made: our data only reveal a clear, usable excess 
at high energies for Revs. 0157, 1093 and 1094 (i.e. 
when HD 5980 is the brightest). This high-energy 
emission was fitted by an absorbed optically-thin 
plasma model (mekal, in Xspec vl 1.2.0): com- 
pared to the spectral properties reported by Naze 
et al. (2002), only changes of the flux level were 
detected; within the confidence intervals, the ab- 
sorption and temperature do not seem to vary in 
a significant way. However, it must be noted that 
the spectra are rather noisy and that only large 
variations would be detected. 

4. Discussion 

The observed variations of the X-ray emission 
from HD 5980 can be qualitatively explained by 
considering the geometry of the CW region in the 
close A+B binary. 

The change in the X-ray flux could be related 
to absorption. Indeed, at phases close to 0.36, the 
opening of the shock cone is directed towards the 
observer and the X-ray emission of the CW region 
is thus seen through the lower density wind of Star 
B (Koenigsberger 2004 and Fig. 1). The lower 
absorption results in an increase of the observed 
X-ray luminosity. At other phases, the very dense, 
slower wind of star A strongly absorbs at least part 
of the CW emission; the width of the lightcurve is 
then directly related to the opening angle of the 
shock cone. Such a scenario has been proposed to 
explain the behaviour of the WR+O binary 7 2 Vel 
(Schild et al. 2004). 

However, the available observations do not yet 
enable us to ascertain that the maximum bright- 
ness really occurs at/near 0=0.36. For example, it 
is conceivable that the count rate actually contin- 
ues to increase beyond this phase, e.g. to peak at 
apastron. A similar variation, proportional to the 
separation between the stars rather than inversely 
proportional, was observed for Cyg OB2 #8A (De 
Becker et al. 2006). In that case, the variation 
was explained by a sudden decrease of the wind 
velocity at periastron because of a strong radia- 
tive inhibition/braking or as a result of perturba- 
tions associated with tidal bulge interactions. For 
HD 5980, such changes of the wind velocity would 
mean that the X-ray emission should not only be 
stronger, but also harder towards apastron, as sug- 
gested by the available data. 


To distinguish between these two scenarios, ad- 
ditional X-ray observations are necessary in order 
to complete the lightcurve near apastron. This 
would enable us to notably determine the posi- 
tion of the brightness peak (apastron or eclipse 
of A by B ?) and the exact shape of this in- 
crease. A detailed spectral analysis is also needed, 
but it requires better spatial resolution than that 
of XMM-Newton in order to cleanly disentangle 
HD 5980 from the superposed SNR, therefore per- 
mitting access to the lower energy range which is 
more affected by the absorption variations. 

In addition, hydrodynamical modelling of the 
fast wind/slow wind interaction region and the 
CW area should be undertaken to quantify their 
respective strength and expected modulation. In 
this context, the lack of strong secular variations 
over 5 years is rather puzzling, since, during this 
interval, the wind density of Star A changed, as 
witnessed in the UV and optical emission lines 
(Koenigsberger 2004). However, the X-ray lu- 
minosity (in the radiative limit) goes as Lx oc 
M x vl © (Stevens et al. 1992), and since the 
eruption the terminal velocity of the wind of star 
A has increased from about 600 km s" 1 to about 
2000 km s^ 1 (Koenigsberger 2004), which off- 
sets the order-of-magnitude decline in the mass 
loss rate over that interval. On the other hand, 
the relative constancy of the X-ray brightness at 
0 = 0.36 despite a significant decline in the mass 
loss from Star A might indicate that Star A’s wind 
has mainly changed in the direction perpendicu- 
lar to the orbital plane than it has in the orbital 
plane, as might be expected if the mass loss from 
star A is not spherically symmetric (Villar-SbafS 
et al. 2003, (i.e. bipolar - see ). Note however 
that the observed P-Cygni UV absorption lines in- 
dicate that at least some of Star A } s wind vari- 
ation occurred near the orbital plane (Koenigs- 
berger 2004, and references therein). 

5. Summary and Conclusions 

An XMM -Newton monitoring campaign of the 
peculiar massive binary HD 5980 has unveiled the 
variations of its X-ray emission. As two observa- 
tions taken 5 years apart (a separation of more 
than 90 binary orbits!) present similar count 
rates, the main source of X-rays cannot be the fast 
wind/slow wind interaction following the LBV-like 
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eruption of Star A, since its emission is expected 
to monotonically and rapidly decrease with time. 
Because individual Wolf-Rayet stars and LBVs 
are only weak X-ray sources, the high-energy ra- 
diation must be associated with the wind-wind 
collision in the close A+B pair, a fact further 
supported by the detection of phase-dependent 
changes. This is the first time that the presence 
of X-ray emitting gas produced by the collision of 
winds in a binary has been confirmed outside our 
Galaxy. 

The X-ray emission of HD 5980 appears mod- 
ulated with phase: a clear increase is observed 
towards the time of the eclipse of Star A by its 
companion. This can be due either to the lower 
absorption inside the shock cone (a situation rem- 
iniscent to that of y 2 Vel) or to a lower wind veloc- 
ity at periastron (similar to the behaviour of Cyg 
OB2 #8 A). Additional data and hydrodynamical 
modelling are now needed to further distinguish 
between these two scenarios. 
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Fig. 1. — Representation of geometry of HD 5980’s 
wind- wind collision region at the 4 orbital phases 
for which XMM -Newton data are available. The 
shape of the shock surfaces was computed under 
the thin shock approximation (Canto et al. 1996) 
and assuming that Ma = lxlO" 4 Mq yr"" 1 , Mb 
= 2xl0” 5 M 0 yr” 1 , v oo (A)=2000 km s^ 1 , and 
v oo (B)=2600 km s -1 . The broken circle around 
star A represents the extent of the wind acceler- 
ating region (Koenigsberger et al. 2006). 



Fig. 2. — Bottom: Evolution of the count rates 
and hardness ratios with phase from XMM- 
Newton observations (see definitions in Table 1). 
Top: Count rate of HD 5980 measured during the 
lOOks Chandra observation, in the same energy 
band. 
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Fig. 3. — Combined MOS images of the area sur- 
rounding HD 5980 in the 0.5-10. keV energy range. 
The data have been binned to get pixels of 2.5” 
and then smoothed with a gaussian of cr=3px. 
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Table 1: Results from the XMM -Newton observ- 
ing campaign. Count rates were evaluated by the 
simple eregionalyse task and are given in the 1.5- 
10. keV band (in units of 10” 3 ctss^ 1 ). Hardness 
Ratios (HR) are defined as with 

M and H the count rates in the 1-2 and 2-10 keV 
bands, respectively. Phases refer to the ephemeris 
of Sterken & Breysacher (1997), while dates are 
Julian Dates minus 2 450000d. 


Rev. 

Date 

Duration 

0 

Count rate MO Si 

Count rate MOS2 

Count rate pn 

HR (MOSl) 

0157 

1835.245 

0.245d 

0.37 

8.26±0.85 

6.95±0.79 

14.92±1.47 

-0.46 ± 0.06 

0357 

2234.643 

0.322d 

0.10 

3.00±0.54 

2.54±0.60 

6.41±1.16 

-0.77 ±0.08 

1093 

3701.863 

0.218d 

0.26 

4.64±0.67 

5.16±0.75 

17.32±1.44 

-0.60 ±0.07 

1094 

3703.807 

0.207d 

0.36 

6.13±0.80 

8.14±0.88 

17.41±1.45 

-0.56 ±0.07 

1100 

3716.125 

0.208d 

1.00 

1.70±0.73 

2.60±0.69 

8.38±1.29 

-0.80 ±0.10 
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